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Magnetization, specific heat, and neutron scattering measurements were performed to study a 
magnetic transition in jarosite, a spin-| kagome lattice antiferromagnet. When a magnetic field 
is applied perpendicular to the kagome plane, magnetizations in the ordered state show a sudden 
increase at a critical field He, indicative of the transition from antiferromagnetic to ferromagnetic 
states. This sudden increase arises as the spins on alternate kagome planes rotate 180° to ferromag- 
netically align the canted moments along the field direction. The canted moment on a single kagome 
plane is a result of the Dzyaloshinskii-Moriya interaction. For H < He, the weak ferromagnetic in- 
terlayer coupling forces the spins to align in such an arrangement that the canted components on 
any two adjacent layers are equal and opposite, yielding a zero net magnetic moment. For H > He, 
the Zeeman energy overcomes the interlayer coupling causing the spins on the alternate layers to 
rotate, aligning the canted moments along the field direction. Neutron scattering measurements 
provide the first direct evidence of this 180° spin rotation at the transition. 

PACS numbers: 75.25.-j, 75.30.Kz, 75.50.Ee, 75.40.Cx 



I. INTRODUCTION 

Frustrated spin systems have received considerable 
attention from both theorists and experimentalists 
due to their unconventional ground states and spin 
dynamics such as spin liquidsfi"— spin nematic;^ and 
spin icei^i Antiferromagnetically coupled spins located 
on the vertices of corner-sharing triangles that make 
up the kagome lattice (Figure [TJa)) are highly sought 
after as being one of the most highly frustrated spin 
systems on a two-dimensional (2D) latticei^i^ Among 
other realizations of the kagome lattice antiferro- 
magnet such as SrCrgGaaOigr^iii volborthiter^"— 
herbertsmithite^-iS and RbaCuaSnFia A— , 
jarosite^^"— is considered one of the most ideal 
systems for the following reasons. First, it consists of 
single and undistorted kagome planes with no structural 
phase transition down to the base temperature. Second, 
a stoichiometrically pure, single-crystal sample can be 
synthesizedi^ 

The general chemical formula of materials in this fam- 
ily is (013)6(804)2, where ^ is a monovalent or bi- 
valent cation such as Na+, K+, Rb+, NHJ, Ag+, T1+, 
HsO"*", or iPb^+, and M is a trivalent cation such as 




FIG. 1: (Color online) (a) The kagome lattice is drawn as 
corner-sharing triangles with the q — ground state spin 
structure. Plus signs inside triangles indicate the sign of the 
vector chirality. (b) The structure of jarosite is shown as 
a ball-and-stick model, highlighting the Fe^^^3(/i-OH)3 bond 
connectivity. 

Fc'^+, Cr''+, or V^^^. For jarosite, the magnetic S — 5/2 
Fe'^"'' ions, which lie at the center of a tilted octahedral 
cage formed by six oxygen atoms, reside at the corners of 
the triangles of the kagome lattice. The kagome planes, 
which are separated by nonmagnetic cations A and sul- 
fate groups, stack up along the c-direction with weak 
interlayer coupling, attesting to the two-dimensionality 
of this system. 
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Even though jarositc has long been regarded as a prime 
model for studying 2D spin frustration, until recently 
the inability to synthesize magnetically pure samples had 
clouded efforts to fully understand the precise magnetic 
behavior of this system. We have developed kinetic con- 
trol of the sample preparation process through oxidation- 
reduction chemistry, which leads to full occupancy of 
the magnetic sites This technique affords mate- 

rials that display consistent magnetic properties. Most 
notably, we observe a prominent antiferromagnetic or- 
dering temperature Tjv narrowly ranging from 56-65 K 
for jarosites with different interplanar cations In ad- 
dition, the nearest-neighbor exchange interaction is con- 
sistent for all jarosites, indicated by similar values of the 
Curie- Weiss temperature of ^ —800 K. This similarity 
in the Curie- Weiss temperature is not surprising, as the 
basic magnetic unit, the Fe^^^3(/x-OH)3 triangle, is struc- 
turally homologous in all jarosites. 

Due to frustration and low dimensionality, ground 
states of a Heisenberg kagomc lattice antiferromagnet 
are infinitely degenerate. Hence, no three-dimensional 
(3D) magnetic order is expected at any non-zero tem- 
perature. However, in jarosite, anisotropic terms and 
small interlayer coupling in the spin Hamiltonian are 
present, and the 3D magnetic long-range order (LRO) 
has been observed at low temperature T < Tj^. Powder 
neutron diffraction measurements on KFe3(OH)6 (804)2 
show that the spins order in a coplanar q = structure 
(Figure ma)) below Tn « 65 K,^^ In the g = struc- 
ture, spins on each triangle form the 120° arrangement 
within the kagome plane, and point either toward (all-in) 
or away from (all-out) the center of the triangle. 

The ordered spins can be decomposed into three sub- 
lattices, with the spins on each triangle labeled Si, S2 
and S3 as shown in Figure [TJa). The vector chirality for 
each triangle can be defined as: 



K,, 



3\/3 



(Si X S2 + S2 X S3 + S3 X Si) 



(1) 



For the coplanar arrangement, this vector is parallel to 
the c-axis with amplitude of -t-1 or -1. A triangle has pos- 
itive (negative) vector chirality if the spins rotate clock- 
wise (counterclockwise) as one traverses around the ver- 
tices of the triangle clockwise. Figure [TJa) shows the 
q = structure with positive vector chirality. Figure [2] 
shows the stacking of the kagome layers along the c- 
direction, where white triangles represent the all-in ar- 
rangement and black triangles represent the all-out ar- 
rangement. The scalar chirality for each triangle can also 
be defined as: 



Ks = Si- (S2 X S3) 



(2) 



The scalar chirality on a single plane is non-zero if spins 
cant out of the plane. We have previously reported the 
field-induced spin-canting transition that corresponds to 
a non-trivial change of the scalar chirality from the net 
Ks = Oto Ks^O state,^ 




H<H 



H>H 



FIG. 2: Spin re-orientations are shown on the second, fourth 
and sixth layers (from the bottom) when H > He- White 
triangles indicate the all-in arrangement and dark triangles 
indicate the all-out arrangement. The shaded area shows the 
volume of the magnetic unit cell. 



The 3D magnetic LRO in jarosite is a result of ferro- 
magnetic interlayer coupling and spin anisotropy arising 
from the Dzyaloshinskii-Moriya (DM) interactionj^i*^i2^ 
The DM interaction, which is a perturbation on the 
Heisenberg spin Hamiltonian, is present if there is no 
inversion center between magnetic sites, and its strength 
is linearly proportional to the spin-orbit couplingi^"— 
It has been experimentally observed in spin frustrated 
perovskite cuprates;^^— the pyrochlore antiferromagnet 
Cu403j^ and molecule-based magnets4^"— The micro- 
scopic calculations of the DM interaction on the kagome 
lattice have been carried out by Elhajal et al^ by apply- 
ing Moriya's formulatioui^ To first approximation, the 
spin Hamiltonian is given by: 



[JiS, -Sj+D,, -S, xS,], 



(3) 



where !]„„ indicates summation over pairs of nearest 
neighbors, Ji is the nearest-neighbor exchange interac- 
tion, and Dij is the DM vector defined in a local axis<^ 
Our spin wave measurements on a single-crystal sample 
provide all relevant spin Hamiltonian parameters for this 
Heisenberg spin modeli^i^i^ 

Figure [3] shows the DM vector, which lies on a mirror 
plane bisecting the bond connecting two magnetic sites, 
with respect to the tilt of the FeOe octahcdra. If one 
octahedron tilts outward, the octahedron that lies right 
on top or below will tilt inward. Therefore, the direction 
of the out-of-planc component of the DM vector changes 
from out-of to into the page or vice versa when one moves 
up or down one layer. On the other hand, the in-plane 
component always points toward the center of the trian- 
gle. Unfortunately, these symmetry considerations alone 
do not determine a precise direction or value of the DM 
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FIG. 3: DM vectors on the kagome lattice. The tilt of the 
octahedra is shown in the lower diagram. The DM vector lies 
within the mirror plane bisecting the bond connecting two 
magnetic sites. In this case, it is assumed that the outward 
(inward) tilt of an octahedron corresponds to the DM vector 
pointing into (out-of) the page. 



vector, which depends on microscopic details, and can 
only be calculated by a microscopic theory. However, 
once the DM vector is calculated for one bond, all other 
DM vectors can be determined using symmetries of the 
crystal. 

The DM interaction can induce a magnetically or- 
dered state at a non-zero temperature by lifting the 
macroscopic ground-state degeneracy, and determines 
the ground-state spin structurei^i^ In jarosite, the tilt 
of the octahedra is about 20° with respect to the crys- 
tallographic c-axis;2L2^ resulting in both out-of-plane 
and in-plane Dp components of the DM vector. Elhajal 
et al. show that the DM interaction can give rise to LRO 
in the kagome lattice, and induce spins to cant out of 
the kagome planes to form an "umbrella" structure. Dz 
confines the spins to lie within the kagome planes, and 
hence effectively takes a role of an easy-plane anisotropy. 
Furthermore, it breaks the symmetry between positive 
and negative vector chirality, leading to the q = spin 
structure. On the other hand. Dp, which takes a role of 
an easy-axis anisotropy, breaks the rotational symmetry 
about the c-axis and gives rise to the all-in-all-out spin 
arrangement. It also forces the spins to cant out of the 
kagome planes to form the umbrella structure of ferro- 
magnetically aligned canted moments giving each layer 
a non-zero magnetic moment i^i^i^ However, in the ab- 
sence of an applied field, the intcrlayer coupling prompts 
the canted moments on neighboring planes to align anti- 
ferromagnetically along the c-axis (Figure [5]) , yielding a 
zero net magnetic moment. 

In this article, we report the results of specific heat. 
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FIG. 4; (Color online) Specific heat of polycrystalline sam- 
ples of deuterated Ag jarosite (red circles), and of non- 
magnetic isostructural compound KGa3(OH)6(S04)2 (blue 
squares) measured at 7/ = T. The KGa3(OH)6(S04)2 data 
are scaled to match the Ag jarosite data at high temperature 
and used to estimate lattice contribution. The inset shows 
the specific heat of Ag jarosite at different magnetic fields. 
Uncertainties where indicated in this article are statistical in 
origin and represent one standard deviation. 



magnetization, and neutron scattering measurements 
on two jarosites, KFe3(OH)6 (804)2 (K jarosite) and 
AgFe3(OH)6(S04)2 (Ag jarosite). In Section II, we dis- 
cuss the specific heat and magnetization measurements, 
and study the field-induced transition between states 
with different spin textures. In Section HI, we investi- 
gate spin re-orientation at the transition using neutron 
scattering. The result provides the first direct evidence 
of the long-expected spin rotation. In Section IV, the 
phase diagrams of Ag and K jarosites are discussed. 



II. SPECIFIC HEAT AND MAGNETIZATION 
MEASUREMENTS 

Details of the synthesis and X-ray crystal structure re- 
finement of Ag jarosite have been previously reported in 
Ref. |2^, whereas those of K jarosite have been reported 
in Refs. [13 and [Sll. The specific heat C of a deuter- 
ated polycrystalline sample of Ag jarosite was measured 
over a 5-150 K temperature range at field strengths vary- 
ing from 0-14 T (Figure |4]). At zero field, a cusp at 
T/v « 59 K indicates the transition to the 3D magnetic 
LRO. When the field is aligned along the c-axis, the cusp 
broadens with increasing field (the inset of Figure U). 
The entropy associated with the magnetic order, which 
is obtained from integrating the magnetic part of C /T 
(after correcting for lattice contribution which was es- 
timated from the non-magnetic isostructural compound 
KGa3(OH)6 (804)2) over the temperature range from 5 K 
to 150 K, represents about 70% of the i?ln6 (where R is 
the molar gas constant) total entropy expected from the 
spin-5/2 system, compared with ~ 50% for K jarositei^ 
This result suggests that short-range correlations have 
already formed at temperature higher than 150 K. 
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FIG. 5: Magnetic susceptibility and specific heat were mea- 
sured on single-crystal K jarosite with magnetic fields up to 
14 T. A field was applied parallel to the c-axis. (a) M/H 
as a function of temperature measured at H —7 T (closed 
squares), 13 T (open squares), and 14 T (open triangles), 
(b) Specific heat as a function of temperature measured at 
H — T (open circles) and 13.7 T (closed circles). The inset 
shows the spin-only specific heat. 



The specific heat of single-crystal K jarosite measured 
with the field of 13.7 T along the c-axis shows two anoma- 
lies, one at 64 K and the other at 50 K (Figure O. The 
former coincides with the transition to the ordered state 
in zero field. The sharp peak is, however, replaced by 
a broad cusp in the presence of the magnetic field. The 
latter occurs at the same temperature where the sharp 
drop in the magnetization is observed, corresponding to 
the spin re-orientation transition between the ferromag- 
netic and antiferromagnetic states. This anomaly is only 
observed in the single-crystal sample, but not in the poly- 
crystalline sample due to powder average and its small 
contribution to the overall specific heat. 

Magnetization measurements on polycrystalline sam- 
ples of non-deuterated and deuterated Ag jarosite were 
performed over a 5-300 K temperature range at field 
strengths varying from 0-14 T. All raw magnetization 
data were corrected for diamagnetism due to core elec- 
trons, and paramagnetic contributions due to free spins 
or impurities. To extract the paramagnetic contribu- 
tions, the low field region of the magnetization M{H) 
for < 9.5 T at T = 5 K was fit to M{H) = 
PiBj{H/T) + P2H + P3, where P„ are empirical pre- 
factors and Bj(x) is the Brillouin function with J = 5/2. 
The average value of Pi upon increasing and decreas- 
ing fields at 5 K is about 0.005 /zb, which indicates that 
about 0.1% of the ¥e^+ spins are non-interacting. This 
result suggests that all magnetic sites are occupied, and 
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FIG. 6: (Color online) (a) Zero field-cool and field-cool (inset) 
magnetic susceptibility of deuterated Ag jarosite. (b) shows 
Tc as a function of field. 



Ag jarosite has virtually no non-interacting spins in the 
ordered state. 

The magnetic susceptibility (x = M / H) measure- 
ments on non-deuterated samples at low field of 100 Oe 
(not shown) show LRO below the Neel ordering tempera- 
ture Tn of 59.6 K.^ At high temperature (T> 150 K), x 
follows the Curie- Weiss law x = C / {T ~ Q'^^^) . A fit to 
this equation yields Q'cw = -800(50) K and C" = 5.4(5) 
cm'^ K/mol Fe. The ordering temperature Tn is greatly 
reduced from the mean field value \&cw\ "^^^ to geomet- 
ric frustration, suggested by a large empirical frustration 
parameter^ / = \Q'cw\IT'n ~ 13. Therefore, the mean 
field theory cannot be applied to determine the effective 
moment /ie//, and the nearest neighbor exchange cou- 
pling Ji- Harris et al^ introduce the corrections to the 
mean field theory to account for geometric frustration us- 
ing the high-temperature series analysis for the kagome 
lattice, which gives = |6cw and C = |C where 

Qcw and C are the mean field Curie- Weiss temperature 
and constant, respectively. Using these corrections, we 
obtain the effective moment /ie// = 6.2(6) fiB and the 
nearest neighbor exchange coupling Ji = 3.9(2) meV. 
For comparison, the values of the nearest neighbor cou- 
pling and the effective moment for K jarosite deduced 
from the Curie-Weiss fit using Harris' corrections are 
Ji = 3.9(2) meV and iieff = 6.3(2) fisr^ respectively. 
We note the similarity of the values of Ji and A*e// for 
K jarosite and Ag jarosite due to the structural homol- 



ogy of the fundamental interacting unit, the Fe^^^3(/i- 
0H)3 triangle. Goodenough-Kanamori rules governing 
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the frontier orbitals of the hydroxide-bridged iron trimer 
predict this antiferromagnetic exchangei ^'''^^ The singly 
occupied dx2_y2 orbitals of high spin Fe^+ overlap with 
the filled sp^ hybrid orbitals of the hydroxide, giving a 
superexchange pathway of cr-synimetry. 

Figure [BJa) shows the magnetic susceptibility of 
deuterated Ag jarosite after subtracting the paramag- 
netic part. All measurements were performed on warm- 
ing at fields between 3 T and 14 T. At high fields, the 
peak at T ss 60 K becomes a broad plateau. The plateau 
is getting broader as the field increases, and extends to 
the base temperature of 5 K for H > 10.5 T. The zero- 
field-cool and field-cool data show the same tempera- 
ture dependence except at low temperature where the 
field-cool data show a rise in the susceptibility, which we 
believe is due to trapped ferromagnetic domains of the 
polycrystalline sample. The susceptibility of the non- 
deuterated sample (not shown) shows the the same tem- 
perature and field dependences but the saturated value 
at high fields (xmax ~ 9 x 10~^ cm^/mol Fe) is about 
10% smaller than that of the deuterated sample. 

Figure [DJa) also shows sharp decreases in the suscep- 
tibility in the ordered state. By taking the derivative 
of the susceptibility with respect to temperature, we 
can define Tc at a peak in the derivatives. Tc, which 
is field-dependent, indicates the transition between the 
ferromagnetic and antiferromagnetic states. Figure IHJb) 
shows Tc as a function of field for deuterated Ag jarosite. 
For T < Tc{H), the system is in the magnetically ordered 
state with no net ferromagnetic moment. In this state, 
the net scalar chirality is zeroi^ For Tc{H) < T < Tn, 
the system is also in the magnetically ordered state, but 
the net scalar chirality is now non-zero and the net fer- 
romagnetic moment is present. From the magnetization 
and specific heat measurements, it is difficult to precisely 
define the ordering temperature Tn at high fields due 
to the lack of well-defined peaks. However, Tjv can be 
readily measured by neutron scattering, which will be 
discussed Section III. 
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FIG. 7: (Color online) (a) Magnetization as a function of field 
at T = 5 K for deuterated Ag jarosite shows hysteresis as the 
field is swept up and down (1 emu = 10~^ A-m^). (b) Tem- 
perature dependences of the critical field He between 5 K 
and 60 K are obtained from magnetization measurements on 
non-deuterated Ag jarosite (blue circles) and deuterated Ag 
jarosite (red squares). Neutron scattering results on deuter- 
ated Ag jarosite measured at BENSC and NCNR are denoted 
by green diamonds and green triangles, respectively. A green 
curve corresponds to a power-law fit to a\T — Tn\^ between 
35 < r < 65 K with exponent /3 = 0.20(1) and Tat = 58.8(3) 
K. Open and closed symbols denote decreasing and increas- 
ing field sweeps, respectively. The inset shows the spin re- 
orientation transition of the first layer from the all-in arrange- 
ment for H < He to the all-out arrangement for H > He- 



A. Critical field, He, and canted moment, AM 

High field magnetization measurements on polycrys- 
talline samples of Ag jarosite were performed up to a 
maximal field of 14 T; the similar measurements on a 
single-crystal sample of K jarosite were previously re- 
ported in Ref. [3l|. Figure Wi^) shows the magnetization 
M at T = 5 K as a function of field after subtracting the 
paramagnetic part. M increases linearly with increasing 
applied field up to 9.75 T, which is followed by a sharp 
increase between 9.75 — 11.5 T. The linear response is 
then recovered for > 11.5 T. Continuing the measure- 
ment with decreasing applied field, M decreases linearly 
to 10.5 T, then drops sharply between 8.75 - 10.5 T. 
The linear behavior then resumes for H < 8.75 T. The 
hysteresis in the bulk magnetization signifies ferromag- 
netic ordering, and is indicative of a first-order tran- 



sition between the antiferromagnetic and ferromagnetic 
states. The anisotropic crystal field of Fe"^"*" ions in a dis- 
torted octahedral geometry yields the hysteresis in the 
M-H curve with a coercive field of approximately 2 T 
at T = 5 K. The critical field He, at which the sudden 
change in magnetization is observed, is defined as a max- 
imum of (dM / dH)\T ■ Upon increasing (decreasing) field. 
He = 10.75 T (8.75 T) at T = 5 K. The same measure- 
ments were conducted at different temperatures to obtain 
Hc{T). The results are shown in Figure [TKb). The line 
corresponds to a fit to the power law HciT) cx jT — T/v|'' 
for a temperature range between 35 K and 65 K giving 
exponent /3 = 0.20(1) and Tn = 58.8(3) K compared 
with /3 = 0.25 and Tn = 64.7 K for K jarosite.-^ 

Additional high-field magnetization measurements on 
a single-crystal sample of K jarosite were performed us- 
ing a nondestructive pulsed magnet, which can gener- 
ate pulsed magnetic fields up to 55 T. at the Interna- 



6 




o 

E 3 



1 - 



(b) 



Li 



10 20 

H [T] 



30 



FIG. 8: (Color online) High-field magnetization measure- 
ments using a pulsed magnet on single-crystal K-jarosite. (a) 
Magnetization curves at several temperatures show a sharp 
rise at He indicated by red arrows and a kink indicated by 
black arrows, (b) dM/dH curves show a clear peak, which 
is used to define He. The data in (a) and (b) are offset for 
clarity. 



tional MegaGauss Science Laboratory, Institute for Solid 
State Physics, The University of Tokyo. These pulsed- 
field measurements are essential to determine the canted 
moment and critical field in K jarosite at low temper- 
ature, since the critical field below 50 K exceeds 14 T, 
the maximal field of our previous measurements Mag- 
netizations were measured by induction using a coaxial 
pick-up coil. A single-crystal sample of K jarosite was 
aligned so that the magnetic field was applied perpendic- 
ular to the kagome plane. The sample was then cooled 
in a liquid *He cryostat down to 4.2 K. Raw data were 
normalized to magnetizations reported in Ref. Isil . 

We also observed the hysteretic behavior, indicative 
of the first-order transition, in K jarosite at low tem- 
perature as shown in the inset of Figure |H1 The coer- 
cive field in K jarosite is, however, much smaller than 
that in Ag jarosite. The hysteresis becomes smaller as 
temperature increases, and is undetectable above 50 K. 
Nevertheless, we cannot rule out the possibility that the 
hysteresis exists up to T/v. In addition to the transition 
at iJc) the magnetization curves show an additional kink 
at Hk > He denoted by black arrows in Figures [H^a) 
and mb). As temperature decreases, Hk and He con- 
verge, and the high-slope region between Hk and He dis- 
appears at the base temperature. This kink was not 
clearly observable in the polycrystalline sample of Ag 
jarosite. It is possible that the prominent first-order tran- 
sition in Ag jarosite suppresses this intermediate phase. 
A similar kink in the M-H curve was also observed in fer- 
romagnetic perovskites, which show a change from first- 
to second-order magnetic phase transitions as a function 
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FIG. 9: (Color online) (a) Magnetization as a function of 
field at different temperatures for deuterated Ag jarosite. The 
H > He range is fit to a linear function with the same slope as 
that of the H < He range. The canted moment AM is defined 
as a t/-intercept. (b) shows temperature-dependence of AM 
for non-deuterated Ag jarosite (blue circles) and deuterated 
Ag jarosite (red squares). AM saturates at 300(30) emu/mol 
Fe (390(30) emu/mol Fe) for the non-deuterated (deuterated) 
sample. Open and closed symbols denote decreasing and in- 
creasing field sweeps, respectively. 



of dopingi^ 

To extract a canted moment AM from the data, we fit 
the high field data {H > He) of the M-H curve to a lin- 
ear function with a fixed slope and extrapolated that line 
to intercept the y-axis to obtain AM (Figure IHl^a) ) . The 
slope was obtained from fitting low field data {H < He) 
to a linear function fixing a constant term to zero. The 
hysteresis in AM persists up to T w 50 K, which is 
in agreement with the hysteresis observed in the criti- 
cal field (Figure EJb)). Below T « 50 K, the canted 
moments measured on increasing field are slightly larger 
than those measured on decreasing field. Figure |9l[b) 
also shows the difference of AM between deuterated and 
non-deuterated samples. This difference could be due to 
a difference in the in-plane component of the DM inter- 
action Dp arising from the replacement of hydrogen by 
deuterium. 

The presence of the canted moment allows us to in- 
vestigate the DM interaction and the intcrlayer coupling 
directly. The sudden change in magnetization is a result 
of the difference in the magnetic moment along the c-axis 
caused by 180° rotation of spins on the alternate layers 
that were previously oppositely canted due to the fer- 
romagnetic interlayer coupling as shown in the diagram 
in Figure [Tjb). At the critical field, the magnetic field 
energy overcomes the ferromagnetic intcrlayer coupling. 



7 



and forces the spins on the alternate layers to rotate 180° 
causing the canted moments to align. 



B. Canting Angle 77, DM Parameters Dp and D^, 
and interlayer coupling J± 

The canting angle, DM parameters, and interlayer cou- 
pling can be calculated from the magnetization results 
{He and AM). AM saturates for temperature below 
30 K, and goes to zero as temperature approaches Tjv. 
The dashed lines in Figure show that AM satu- 
rates at AM(0) 300(30) emu/mol Fc = 0.054(5) hb 
(390(30) emu/mol Fc = 0.070(5) ^s) for non-dcuterated 
(deuterated) Ag jarosite. However, since the data were 
measured on polycrystalline samples, AM is multiplied 
by three to take into account the powder average. Given 
the expected moment M = 5 fj-B for S = 5/2, the 
canting angle at zero temperature rj is defined as rj = 

sin-i (^1^) , which gives 77 = 1.8(2)° (jy = 2.4(2)°) for 

the non-deuterated (deuterated) sample. For K jarosite, 
AM(0) = 650(80) emu/mol Fe and hence tj = 1.3(2)°, 
which is lower than the value of 1.9(2)° obtained from 
the spin wave measurements 1^ This difference could be 
due to thermal or quantum fluctuations, which reduce 
the staggered moment. The similar measurements of the 
spin canting in La2Cu04, a parent compound of high- 
Tc superconductors, yields a much smaller value of the 
canting angle (77 = 0.17°)i^ 

Since the spin canting is a result of the DM interaction, 
one can write down the expression for rj in terms of the 
DM parameters. Elhajal et al. derived the following 
relation between 77 and the DM parameters Dp and ^ 



1] 



1 



tan 



-2Dr. 



(4) 



.V3Ji-Dz 

We note that signs of Dp and D^ are reversed and dif- 
ferent from Ref. [Sy to be consistent with the DM vector 
defined in this article and in Ref. IH. The absolute- value 
indicates that the sign of the canting angle with respect 
to the in-plane order cannot be uniquely determined from 
the DM vector alone. In fact, for a given value of Dp with 
Dz < 0, there exist two lowest-energy degenerate states 
with opposite canting and different in-plane orders cor- 
responding to the Kg > and Kg < states At zero 
field, this chiral symmetry is spontaneously broken at Tjv, 
and the interlayer coupling selects the spin structure on 
neighboring planes, giving rise to the 3D magnetic LRO. 
Equivalently, applying a magnetic field along the c-axis 
will also lift this degeneracy. 

Yildirim and Harris have also derived a similar ex- 
pression for the canting angle induced by the DM 
interaction)^ 



1 . 

— sm 
2 



-2D„ 



where J = Ji -I- J2 . Both Eq. |4] and Eq. [5] approach the 
same limiting value for a small canting angle. In order to 



determine both in-plane Dp and out-of plane Dz compo- 
nents of the DM interaction, we also need the value of the 
energy gap of the lifted zero energy mode, and therefore 
rely on previous single-crystal neutron scattering mea- 
surements on K jarosite<^ The spin wave calculations by 
Yildirim and Harris give the following expression for the 
two energy gaps around 7 meV at F pointi^ 



hujQ = s 



iDl 



\%Di - eVsjDz ± 



(6) 



The energy splitting of these two gaps is about 0.1 meV, 
where the energy of the lifted zero energy mode has a 
smaller value. This value is independent of the sign of 
Dp due to the ± sign in front of the last term. \Dp\ and 
Dz can be calculated by solving Eqs.[S] and [5] providing 77 
and hajQ. Given hujQ = 6.7(1) meV from the single-crystal 
neutron scattering measurements^ and using the data 
shown in Figure IIOI we estimated the energy gap for Ag 
jarosite to be 6.2(4) meV. Therefore, after solving Eqs.[5] 
andm we obtain \Dp\ = 0.18(3) meV (0.24(3) meV), and 
Dz = -0.17(3) meV (-0.17(3) meV) for non-deuterated 
(deuterated) Ag jarosite (see Table H]). Here, we assume 
that J is the same for Ag jarosite and K jarosite as evi- 
denced by the the same Curie- Weiss temperature. 

In addition, the knowledge of AM(0) and -ffc(O) from 
the magnetization measurements allows us to calculate 
the interlayer coupling J±. Similar to AM, He ap- 
pears to saturate at low temperature to the average 
value of Hc{0) = 9.5(1.0) T. As temperature increases. 
He monotonically decreases and goes to zero as tem- 
perature approaches T^. The field-induced transition 
is a result of the competition between the ferromag- 
netic interlayer coupling and the Zeeman energy. Us- 
ing the equation Jj_ = -AM(0) • He{0)/S^^^ we ob- 
tain J_L = -0.014(3) meV (-0.018(4) meV) for the non- 
deuterated (deuterated) Ag jarosite. The minus sign in- 
dicates that the interlayer coupling is ferromagnetic. For 
K jarosite, using the canting angle of 1.9(2)° that gives 
AM(0) = 0.17(2)/iB from the spin wave measurements, 
or AM(0) = 0.12(2)/iB from the pulsed-filed magnetiza- 
tion measurements, we obtain J± = —0.027(4) meV or 
-0.019(4) meV, respectively (see Table The magni- 
tude of J± is hundreds times smaller than the magnitude 
of Ji , which attests to the two dimensionality of this spin 
system. 



III. NEUTRON SCATTERING 
MEASUREMENTS 

Inelastic neutron scattering measurements in zero field 
were performed using the thermal triple-axis spectrom- 
eters BT7 and BT9 at the NIST Center for Neutron 
Research (NCNR) with the final energy fixed at 14.7 
meV to measure spin wave excitations on deuterated 
polycrystalline Ag and K jarosites, respectively. The 
horizontal coUimations of open— 50'— S— 60' — 120' and 
40' - 56'-S-52'-open were employed at BT7 and BT9, 



8 



TABLE I: Critical fields, canted moments at zero temperature, canting angles, and spin Hamiltonian parameters for K jarosite 
and Ag jarosite calculated from magnetization and neutron scattering data. 



Compound 


Hc(0) (T) 


AM(0) (ms) 




\Dp\ (meV) 


(mcV) 


Jc (mcV) 


hojo (meV) 


AgFe3(OX)6(S04)2 

I. X = H 

II. X = D 


9.5(1.0) 
9.5(1.0) 


0.16(2) 
0.21(2) 


1.8(2) 
2.4(2) 


0.18(3) 
0.24(3) 


-0.17(3) 
-0.17(3) 


-0.014(3) 
-0.018(4) 


6.2(4) 
6.2(4) 


KFe3(OH)6(S04)2 

I. Neutron scattering^ 

II. Magnetization 


16.8(5) 


0.12(2) 


1.9(2) 
1.3(2) 


0.197(2) 


-0.196(4) 


-0.027(4) 
-0.019(4) 


6.7(1) 



25 



^ 20 
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FIG. 10: (Color online) Inelastic neutron scattering measure- 
ments were performed using the triple-axis spectrometers at 
NCNR. The constant-Q scan at \Q\ = 1.7 A on deuterated Ag 
jarosite (red squares) was measured at BT7, and the constant- 
Q scan at |Q| = 1.5 A on deuterated K jarosite (blue cir- 
cles) was measured at BT9. The scans show powder-average 
spin-wave excitations at energy transfer of 7.9(4) meV and 
8.5(3) meV for Ag jarosite and K jarosite, respectively. The 
K jarosite data were normalized so that they lie on top of the 
Ag jarosite data. Solid lines show the fits to the spin wave 
excitations described in Ref. l34l . convoluted with the instru- 
mental resolution function. 
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respectively. Pyrolytic graphite (PG) crystals were uti- 
lized to monochromate and analyze the incident and scat- 
tered neutron beams. PG filters were placed in the scat- 
tered beam to reduce higher-order contamination. The 
samples were cooled using a closed cycle ^He cryostat to 
T = 5.3 K at BT7 and to T = 14 K at BT9. The energy 
scans in Figure [10] at \Q\ = 1.7 A for Ag jarosite (red 
squares) and \Q\ = 1.5 A for K jarosite (blue circles) 
show powder-average spin wave excitations of the lifted 
zero-energy mode^ which is weakly dispersive, at energy 
transfer of 7.9(4) meV and 8.5(3) mcV, respectively. 

To further study the magnetic transition at i?c, clas- 
tic neutron scattering measurements on a polycrystalline 
sample in a high magnetic field were performed using the 
triple-axis spectrometer El, the diffractometer E4 at the 
Berlin Neutron Scattering Center (BENSC), Helmhohz 
Centre Berlin for Materials and Energy (formerly Hahn- 
Meiter-Institut Berlin), with the incident neutron wave- 
length of 2.43 A {huj = 13.9 meV), and the triple- 



FIG. 11: (Color online) Elastic neutron scattering measure- 
ments on a polycrystalline sample of deuterated Ag jarosite 
were performed using the triple-axis spectrometer El. (a) \Q\- 
scans at T = 3 K show two magnetic Bragg peaks in zero field 
(red circles) at (1, 1, 3/2) and (0, 1, 7/2), and two magnetic 
peaks at = 14 T (blue squares) at (1, 1, 0) and (0, 1, 7/2). 
The structural (1, 0, 4) nuclear peak is denoted by the green 
arrow, (b) The field dependence of the (1, 1, 0) peak intensity 
shows the hysteretic phase transition to the net Ks ^ state 
at high field. 



axis spectrometer BT7 at NCNR with the incident neu- 
tron energy of 14.7 meV. The horizontal collimations of 
40'-80'-S-40'-40' and 40'-opcn-opcn were employed 
at El and E4, respectively, and the horizontal collima- 
tions of open-50'-S-60' - 120' were employed at BT7. 
PG crystals were utilized to monochromate the incident 
neutron beam at El, E4, and BT7, and to analyze the 
scattered neutron beam at El and BT7. PG filters were 
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placed in the incident beam at El and E4, and two 
PG filters were placed before and after the sample at 
BT7 to reduce higher-order contamination. At BENSC, 
deuterated polycrystalline Ag jarosite of mass 5.5 g was 
mounted inside the VM-1 magnet with the highest field 
of 14.5 T, and at NCNR 5.4 g of the polycrystalline sam- 
ple was mounted inside the 11 T magnet. The sample 
was then cooled using a ''He cryostat. 

At zero field, the canted moments on successive kagome 
planes are equal and opposite. It takes six layers to form 
a magnetic unit cell instead of three layers for a struc- 
tural unit cell (Figure [5]). This causes the doubling of a 
magnetic unit cell along the c-axis, resulting in magnetic 
Bragg peaks at a half integer value along the 001 direc- 
tion. For H > He, the spins on alternate layers rotate 
180° changing the spin arrangement from all-in to all-out 
and vice versa. Figure [2] shows such spin re-orientations 
on the second, fourth and sixth kagome planes. In this 
state, the magnetic unit cell contains only three kagome 
planes, and is the same size as the structural unit cell, 
resulting in the magnetic Bragg peaks at a whole integer 
value along the OOZ direction. On the other hand, if the 
transition to the ferromagnetic state is caused by simply 
flipping the canted components on the alternate layers, 
the positions of the magnetic Bragg peaks will remain at 
the half integer value. 

In Figure [Til |Q|-scans at T = 3 K measured on a 
deuterated polycrystalline sample of Ag jarosite at El re- 
veal the shift of the magnetic Bragg peaks expected from 
the 180° spin rotation. To overcome the intcrlayer cou- 
pling and induce the 180° spin rotation, an applied field 
greater than the critical field is required along the crys- 
tallographic c-axis. At zero field, we observe two mag- 
netic Bragg peaks at (0, 1, 7/2) and (1, 1, 3/2) in this 
IQI range. At H = U T {H > He), the peak at (1, 1, 
3/2) disappears, while a new magnetic peak at (1, 1, 0) 
emerges. The intensity at (0, 1, 7/2) does not change 
because the magnitude of the magnetic field component 
along the c-axis is less than the critical field. Since the 
applied field is perpendicular to the scattering plane, at 
(1, 1, 0) the c-axis is along the vertical direction, hence 
along the direction of the applied field. On the other 
hand, at (1, 1, 3/2) and (0, 1, 7/2) the angles between 
the field direction and the c-axis are 18° and 53°, respec- 
tively. Therefore, the field component along the c-axis 
is reduced by about 5% at (1, 1, 3/2) and about 40% 
at (0, 1, 7/2). With the applied field of 14 T, the field 
strengths along the c-axis at (1, 1, 3/2) and (0, 1, 7/2) 
are 13.3 T and 8.4 T, respectively, while the critical field 
atT = 3Kisl0.5T. Thus, we observe a decrease in 
intensity at (1, 1, 3/2), but no change in intensity at 
(0, 1, 7/2) (Figure [TlT a)). The critical field was mea- 
sured using the data taken at (1, 1, 0) where the applied 
field is perfectly along the c-axis, and the change in the 
intensity is largest. The field dependence of the intensity 
at (1, 1, 0) (Figure [TTJb)) illustrates the hysteresis con- 
sistent with the magnetization data. Furthermore, the 
critical fields obtained from the neutron scattering mea- 
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FIG. 12: (Color online) (a) cj-scans around (1, 1, 0) at three 
different magnetic fields at T = 54 K measured at El on 
single-crystal K jarosite. The magnetic field was applied along 
the c-axis, which is normal to the scattering plane. The result 
shows the change in the integrated intensity of the (1, 1, 0) 
refiection as the field exceeds the critical field He- (b) The 
integrated intensity as a function of field measured at T = 54 
K shows a sudden change in intensity at the critical field. 
The inset shows the derivative of the integrated intensity as 
a function of field. 



surements are in good agreement with the magnetization 
results as shown in Figure [Tljb). 

To investigate this magnetic transition on a single- 
crystal sample, the elastic neutron scattering measure- 
ments on single-crystal K jarosite were performed using 
the triple-axis spectrometers El at BENSC and BT7 at 
NCNR with the same instrumental set-up as described 
above. Additional neutron diffraction was performed at 
the diffractometer E4 at BENSC with the incident neu- 
tron wavelength of 2.45 A using a position sensitive de- 
tector. At BENSC, the single-crystal sample (mass « 
20 mg) was mounted inside the VM-1 magnet, and at 
NCNR the same sample was mounted inside the 11 T 
magnet. The sample was cooled using a ^He cryostat. 
It was aligned in the (hkO) plane with the magnetic field 
applied perpendicular to the scattering plane. The inten- 
sity of the (1, 1, 0) refiection was measured as a function 
of field and temperature. At NCNR, we measured the 
intensities of both (1, 1, 0) and (1, 1, 3/2) reflections to 
construct a complete phase diagram. To measure the (1, 

1, 3/2) reflection, the sample was re-aligned on the (2, 

2, 3) structural Bragg peak in air, and the intensity of 
the (1, 1, 3/2) reflection was measured inside the 11 T 
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FIG. 14: (Color online) K jarosite phase diagram as a func- 
tion of field and temperature. Open circles show the results 
of the magnetization measurements. Yellow down triangles 
show the results of the neutron scattering measurements at 
BENSC. Red up triangles and purple squares show the re- 
sults of the neutron scattering measurements at NCNR. The 
critical fields obtained from the pulsed magnet at ISSP are 
denoted by closed and open diamonds for increasing and de- 
creasing field sweeps, respectively. Blue circles denote fields 
at which a kink in magnetization was observed as shown by 
black arrows in Figure [S] A contour plot shows the intensity 
of the (1, 1, 0) magnetic Bragg peak measured at E4. A solid 
line is a power law fit with exponent /? = 0.25. Dotted lines, 
which serve as guides to the eye, denote the phase lines. 



T[K] 



FIG. 13: (Color online) Neutron scattering intensities as a 
function of temperature were measured on a single- crystal 
sample of K jarosite (a) at Q =(1, 1, 0) and H — 11 T 
and (b) at Q =(1, 1, 3/2) and H = 6.7 T. A dotted arrow 
denotes a crossover phase line between the 3D ordered state 
with Ka 7^ and the 2D short-range correlated, disordered 
state. Dotted lines are guides to the eye. (c) shows the field 
dependence of the (1, 1, 0) magnetic Bragg intensity. 



magnet at low temperature. As previously discussed, the 
magnetic field is along the c-axis for (1,1,0) and about 
18° with respect to the c-axis for (1, 1, 3/2). The val- 
ues of magnetic fields in Figures and [T2] indicate the 
component of the field along the c-axis. 

Figure I12f a) shows w-scans measured on a single- 
crystal sample at El around the (1, 1, 0) reflection at 
three different fields at T = 54 K. The data are fit to a 
Gaussian. Figure IT^ b) shows integrated intensity, which 
is extracted from the fit in Figure I12f a) , as a function 
of field at T = 54 K. The inset shows the derivative of 
the intensity, where the critical field is defined at the 
maximum. The same measurements were performed at 
several temperatures to obtain Hc{T). Yellow down tri- 
angles in Figure [T3] show the critical field as a function 



of temperature, which is consistent with the magnetiza- 
tion measurements shown by open circles and diamonds. 
Figures I13f a) and I13f b) show neutron scattering inten- 
sities as a function of temperature and field measured 
at NCNR. The magnetic scattering intensity at (1, 1, 0) 
(Figure [Of a)) is non-zero in the ferromagnetic state, and 
goes to zero above T/v, where spins are disordered, and 
below Tc, where spins are in the antiferromagnetic state. 
Therefore, by measuring the intensity at (1, 1, 0), one 
can extract both Tjy and Tc- The intensity at (1, 1, 3/2) 
(Figure [TM b)). on the other hand, is non-zero in the an- 
tiferromagnetic state, and goes to zero above T/v or Tc- 
The arrows in Figures [TWa) and fOf b) indicate Tc and 
T/v from the (1, 1, 0) data and the (1, 1, 3/2) data. In 
Figure [Til Tc{H) obtained from the (1, 1, 0) data is de- 
noted by red up triangles, while that obtained from the 
(1, 1, 3/2) data is denoted by purple squares. In addi- 
tion, the field dependence of the (1, 1, 0) scattering in- 
tensity measured at E4 (Figure IT^ c) and a contour plot 
in Figure fT4)) shows sharp drops on the low temperature 
side consistent with the spin re-orientation transition ob- 
served in the magnetization measurements. On the high 
temperature side, a crossover between the 3D long-range 
ordered and 2D short-range correlated phases was ob- 
served. 
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FIG. 15: (Color online) Ag jarosite phase diagram as a func- 
tion of field and temperature. Red circles, denoting the first- 
order phase transition, show the average values of the critical 
fields from the magnetization measurements. Blue triangles, 
denoting the second-order phase transition, show the results 
of the magnetization measurements. Open circles show the 
NCNR neutron scattering results. A solid line is a power 
law fit with exponent = 0.20. A dotted line shows an ex- 
pected phase boundary between the net Ks ^ state and 
disordered state. A green arrow indicates the possibility of a 
multi-critical point. 



IV. PHASE DIAGRAMS 

For Ag jarosite, the results of magnetization and neu- 
tron scattering measurements are summarized in the H- 
T phase diagram shown in Figure In the Kg = Q 
phase, the canted moments on adjacent kagome planes 
point in the opposite directions due to the weak inter- 
laycr coupling, resulting in a zero net magnetic moment. 
However, at the critical field iJc, the Zeeman energy over- 
comes the interlayer coupling, and the spins on the alter- 
nate planes rotate 180°, resulting in the reversal of the 
canted moments. In the magnetization measurements, 
the transition results in the abrupt increase of the mag- 
netization at the critical field. In the neutron scattering 
measurements, this transition shifts the magnetic Bragg 
peaks. 

At low temperature, both magnetization and neutron 
scattering intensity show the hysteretic behavior as a 
function of field. The hysteresis, which is indicative of the 
first-order transition, exists for T < 50 K. The magnetic 
transition line changes to second order for T ^ 50 K, rais- 
ing a possibility of the multi-critical point at T w 50 K. 
The hysteresis was also observed in K jarosite with a 
much smaller coercive field, which becomes too small to 
detect above 50 K. However, we note that we cannot rule 
out the possibility that this first-order transition phase 
line could extend up to Tjy. 

The phase diagram of K jarosite is illustrated in Fig- 
ure 1141 The neutron scattering measurements show the 
field dependence of the 3D magnetic ordering tempera- 
ture T/v. At zero field, the 3D magnetic LRO is a result 



of the combined effect of the DM interaction and inter- 
layer coupling. The DM interaction induces the 2D long 
length-scale correlations in the planes, and the interlayer 
coupling propagates spin correlations in the third direc- 
tion. In the presence of a magnetic field, since the field 
favors the net Ks 7^ state while the interlayer coupling 
favors the net Ks = state, they will compete with each 
other. In the low field region, this competition leads to 
a decrease of T^r. In the high field region, on the other 
hand, the Zeeman energy overwhelms the interlayer cou- 
pling, and takes over its role, resulting in an increase of 
Tjv as the field increases. In Figure [T31 the contour plot 
of the (1, 1, 0) magnetic Bragg intensity measured at E4 
shows a sharp transition at {Tc,Hc), which corresponds 
to the spin re-orientation transition, and an indistinct 
phase line between the 3D ordered state with Ks ^ 
and the 2D short-range correlated, disordered state. At 
_ff « 5 T, the two phase lines meet, which is a special 
point where the effects of the interlayer coupling and 
magnetic field balance against each other, and the 3D 
magnetic LRO is solely induced by the 2D correlations in 
the planes. The lack of a well-defined peak in the specific 
heat at a finite field suggests that the phase boundary 
between the 3D K ^ ordered state and the 2D corre- 
lated, disordered state is a crossover rather than a true 
phase transition. The magnetic field applied along the 
c-axis breaks the scalar chirality symmetry, prohibiting 
further spontaneous breaking of the spin symmetry. This 
scalar chirality symmetry corresponds to the two oppo- 
sitely canting degenerate states. The spontaneous break- 
ing of this Z2 symmetry, which corresponds to the phase 
transition from the 2D short-range correlation state and 
the 3D long-rang magnetic ordered state, can only occur 
at zero field. 



V. CONCLUSIONS 

We have investigated the spin canting due to the DM 
interaction in jarosite using specific heat, magnetization, 
and neutron scattering measurements. From the ob- 
served critical field and canted moment obtained from 
our magnetization measurements, we were able to calcu- 
late the canting angle rj, the DM parameters \Dy\ and 
Dz, and the interlayer coupling constant J±. The in- 
plane component of the DM interaction forces the spins 
on a triangle to cant out of the kagome planes to form 
the umbrella structure of ferromagnetically aligned mo- 
ments within the layers. However, at zero field, the weak 
ferromagnetic interlayer coupling prompts the spins on 
neighboring planes to align in such an arrangement that 
the canted moments are equal and opposite. At the crit- 
ical field, the Zeeman energy overcomes the interlayer 
coupling causing the spins on the alternate layers to ro- 
tate 180°. This 180° spin rotation reverses the canted 
moments on the alternate layers causing the sudden in- 
crease in the magnetization. This field-induced spin 
re-orientation transition corresponds to the non-trivial 
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change in the spin-texture. In particular, the transition 
yields a net, non-zero value for the scalar chirality. Our 
elastic neutron scattering measurements in a high mag- 
netic field provide the first direct evidence of the 180° 
spin rotation at the transition. 



ported by the NSF under Grant No. DMR 0239377, and 
in part by the MRSEC program under Grant No. DMR 
02-13282. 



Acknowledgments 

The authors would like to thank P. Ghaemi and T. 
Senthil for many fruitful discussions. This work was sup- 



* Electronic address: lsckittiwit@ma hidol.ac.tIil 
t Electronic address: younglee@mit.edul 

^ P. W. Anderson, Mater. Res. Bull 8, 153 (1973). 
^ T. Koretsune and M. Ogata, Phys. Rev. Lett. 89, 116401 
(2002). 

^ J. T. Chalker and J. F. G. Eastmond, Phys. Rev. B 46, 
14201 (1992). 

* P. Lecheminant, B. Bernu, C. Lhuillier, L. Pierre, and 
P. Sindzingre, Phys. Rev. B 56, 2521 (1997). 

^ J. T. Chalker, P. C. W. Holdsworth, and E. F. Shender, 

Phys. Rev. Lett. 68, 855 (1992). 
® A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan, 

and B. S. Shastry Nature 399, 333 (1999). 

J. Snyder, J. S. Slusky, R. J. Cava, and P. Schiffer, Nature 

413, 48 (2001). 
^ A. P. Ramirez, Annu. Rev. Mater. Sci. 24, 453 (1994). 
^ J. E. Greendan, J. Mater. Chem. 11, 27 (2001). 

A. P. Ramirez, G. P. Espinosa, and A. S. Cooper, Phys. 

Rev. Lett. 64, 2070 (1990). 

C. Broholm, G. Aeppli, G. P. Espinosa, and A. S. Cooper, 
Phys. Rev. Lett. 65, 3173 (1990). 

Z. Hiroi, M. Hanawa, N. Kobayashi, M. Nohara, H. Takagi, 
Y. Kato, and M. Takigawa, J. Phys. Soc. Jap. 70, 3377 
(2001). 

M. Yoshida, M. Takigawa, H. Yoshida, Y. Okamoto, and 

Z. Hiroi, Phys. Rev. Lett. 103, 77207 (2009). 
" F. Bert, D. Bono, P. Mendels, F. Ladieu, F. Due, J. C. 

Trombe, and P. Millet, Phys. Rev. Lett. 95, 087203 (2005). 

M. Shores, E. Nytko, B. Bartlett, and D. Nocera, Journal 

of the American Chemical Society 127, 13462 (2005). 

J. S. Helton, K. Matan, M. P. Shores, E. A. Nytko, B. M. 

Bartlett, Y. Yoshida, Y. Takano, A. Suslov, Y. Qiu, J.-H. 

Chung, et al.. Physical Review Letters 98, 107204 (2007). 
" J. S. Helton, K. Matan, M. P. Shores, E. A. Nytko, B. M. 

Bartlett, Y. Qiu, D. G. Nocera, and Y. S. Lee, Phys. Rev. 

Lett. 104, 147201 (2010). 

P. Mendels, F. Bert, M. A. de Vries, A. Olariu, A. Harrison, 

F. Due, J. C. Trombe, J. Lord, A. Amato, and C. Baines, 

Phys. Rev. Lett. 98, 077204 (2007). 
1^ S.-H. Lee, H. Kikuchi, Y. Qiu, B. Lake, Q. Huang, 

K. Habicht, and K. Kiefer, Nature Materials 6, 853 (2007). 
2" K. Morita, M. Yano, T. Ono, H. TANAKA, K. Fujii, 

H. Uekusa, Y. Narumi, and K. Kindo, J. Phys. Soc. Jpn. 

77, 043707 (2008). 

T. Ono, K. Morita, M. Yano, H. Tanaka, K. Fujii, 
H. Uekusa, Y. Narumi, and K. Kindo, Phys. Rev. B 79, 
174407 (2009). 



K. Matan, T. Ono, Y. Fukumoto, T. J. Sato, J. Yamaura, 
M. Yano, K. Morita, and H. Tanaka, Nature Physics 
(2010). 

A. S. Wills, Phys. Rev. B 63, 064430 (2001). 

T. Inami, M. Nishiyama, S. Maegawa, and Y. Oka, Phys. 

Rev. B 61, 12181 (2000). 
■^^ M. Nishiyama, S. Maegawa, T. Inami, and Y. Oka, Phys. 

Rev. B 67, 224435 (2003). 
2^ S.-H. Lee, C. Broholm, M. F. Collins, L. Heller, A. P. 

Ramirez, C. Kloc, E. Bucher, R. W. Erwin, and N. Lace- 

vic, Phys. Rev. B 56, 8091 (1997). 

D. Grohol, D. G. Nocera, and D. Papoutsakis, Phys. Rev. 
B 67, 064401 (2003). 

B. M. Bartlett and D. G. Nocera, J. Am. Chem. Soc. 127, 
8985 (2005). 

23 F. C. Coomer, A. Harrison, G. S. Oakley, J. Kulda, J. R. 

Stewart, J. A. Stride, B. FAK, J. W. Taylor, and D. Visser, 

J Phys-Condens Mat 18, 8847 (2006). 
^° W. G. Bisson and A. S. Wills, J. Phys.-Condens. Mat. 20, 

452204 (2008). 

D. Grohol, K. Matan, J. H. Cho, S. H. Lee, J. W. Lynn, 
D. G. Nocera, and Y. S. Lee, Nature Materials 4, 323 

(2005) . 

^2 D. G. Nocera, B. M. Bartlett, D. Grohol, D. Papoutsakis, 
and M. P. Shores, Chem. Eur. J. 10, 3850 (2004). 
A. S. Wills, Phys. Rev. B 63, 064430 (2001). 
K. Matan, D. Grohol, D. G. Nocera, T. Yildirim, A. B. 
Harris, S. H. Lee, S. E. Nagler, and Y. S. Lee, Phys. Rev. 
Lett. 96, 247201 (2006). 

T. Yildirim and A. B. Harris, Phys. Rev. B 73, 214446 

(2006) . 

M. Elhajal, B. Canals, and C. Lacroix, Phys. Rev. B 66, 
014422 (2002). 

I. E. Dzyaloshinskii, J. Phys. Chem. Solids 4, 241 (1958). 
T. Moriya, Phys. Rev. 120, 91 (1960). 
D. Coffey, T. M. Rice, and F. C. Zhang, Phys. Rev. B 44, 
10112 (1991). 

L. Shekhtman, O. Entin-Wohlman, and A. Aharony, Phys. 
Rev. Lett. 69, 836 (1992). 

A. Aharony, O. Entin-Wohlman, and A. B. Harris, NATO 

ASI Series, Series E: Applied Science 349, 281 (1998). 
^■^ O. Cepas, K. Kakurai, L. P. Regnault, T. Zinian, J. P. 

Boucher, N. Aso, M. Nishi, H. Kageyama, and Y. Ueda, 

Phys. Rev. Lett. 87, 167205 (2001). 
'''^ O. Cepas, T. Sakai, and T. Ziman, Prog. Theor. Phys. 

Suppl. 145, 43 (2002). 

L. B. Steren, M. Tovar, and S. B. Oseroff, Phys. Rev. B 



13 



46, 2874 (1992). 

M. Hucker, V. Kataev, J. Pommer, U. Ammerahl, 
A. Revcolevschi, J. M. Tranquada, and B. Buchner, Phys. 
Rev. B 70, 214515 (2004). 

R. K. Li, R. Kremer, and J. Maier, J. Solid State Chem. 
146, 488 (1999). 

M. Elhajal, B. Canals, and C. Lacroix, J. Phys.: Condens. 
Matter 16, S917 (2004). 

M. A. Girtu, C. M. Wynn, W. Fujita, K. Awaga, and A. J. 

Epstein, Phys. Rev. B 61, 4117 (2000). 

R. Blinc, K. Pokhodnia, P. Cevc, D. Arcon, A. Omerzu, 

D. Mihailovic, P. Venturini, L. Golic, Z. Trontelj, 

J. Luznik, et al., Phys. Rev. Lett. 76, 523 (1996). 

Z. Wang, B. Zhang, T. Otsuka, K. Inoue, H. Kobayashi, 

and M. Kurmoo, J. Chem. Soc, Dalton Trans, p. 2209 

(2004). 

J. M. Rawson, A. Alberola, H. El-Mkami, and G. M. 

Smith, J. Phys. Chem. Solids 65, 727 (2004). 

P. A. Joy and S. Vasudevan, J. Chem. Phys 99, 4411 



(1993). 

K. Matan, J. S. Helton, D. Grohol, D. G. Nocera, S. Waki- 
moto, K. Kakurai, and Y. S. Lee, Physica B 404, 2529 
(2009). 

^'^ M. Elhajal, B. Canals, and C. Lacroix, Physica B 312, 716 
(2002). 

^■"^ R. Ballon, B. Canals, M. Elhajal, L. C, and A. Wills, J. 

Magn. Magn. Mater. 262, 465 (2003). 
'^^ A. B. Harris, C. Kallin, and A. J. Berlinsky, Phys. Rev. B 

45, 2899 (1992). 
" J. B. Goodenough, J. Phys. Chem. Solids 6, 287 (1958). 
'^^ J. Kanamori, J. Phys. Chem. Solids 10, 87 (1959). 

J. Mira, J. Rivas, F. RivaduUa, C. Vazquez- Vazquez, and 

M. A. Lopez-Quintela, Phys. Rev. B 60, 2998 (1999). 
^° T. Thio, T. R. Thurston, N. W. Preyer, P. J. Picone, M. A. 

Kastner, H. P. Jenssen, D. R. Gabbe, C. Y. Chen, R. J. 

Birgeneau, and A. Aharony, Phys. Rev. B 38, 905 (1988). 



